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A nascent neutron star resulting from stellar collapse is a prodigious source of neutrinos of all flavors. While the 
most basic features of this neutrino emission can be estimated from simple considerations, the detailed simulation 
of the neutrinos' decoupling from the hot neutron star is not yet computationally tractable in its full glory, 
being a time-dependent six-dimensional transport problem. Nevertheless, supernova neutrino fluxes are of great 
interest in connection with the core-collapse supernova explosion mechanism and supernova nucleosynthesis, and 
as a potential probe of the supernova environment and of some of the neutrino mixing parameters that remain 
unknown; hence a variety of approximate transport schemes have been used to obtain results with reduced 
dimensionality. However, none of these approximate schemes have addressed a recent challenge to the conventional 
wisdom that neutrino flavor mixing cannot impact the explosion mechanism or r-process nucleosynthesis. 



1. CORE-COLLAPSE SUPERNOVAE 

Stellar collapse results from the evolution of 
a massive star. How far a star can burn de- 
pends on its total mass M; but because of cool- 
ing by neutrino emission from carbon burning on- 
wards, the core of a massive star eventually be- 
comes degenerate, be it composed of O/Ne/Mg 
(8 M Q < M < 10 M Q ) or Fe (M > 10 M ). 
The core becomes unstable as the ashes of nu- 
clear burning pile up and it approaches the Chan- 
drasekhar mass. Its inner portion undergoes ho- 
mologous collapse (velocity proportional to ra- 
dius), and the outer portion collapses superson- 
ically. Electron capture on nuclei is one proxi- 
mate cause of instability, producing v e that es- 
cape freely until densities in the collapsing core 
become so high that even neutrinos are trapped. 

Collapse halts with a 'bounce' soon after the 
inner core exceeds nuclear density, when a shock 
wave forms at the boundary between the homol- 
ogous and supersonically collapsing regions. The 
shock begins to move out, but stalls at some dis- 
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tance beyond the surface of the newly-born neu- 
tron star as neutrino emission and endothermic 
dissociation of heavy nuclei sap energy from the 
shocked material. 

Neutrinos dominate the energetics of post- 
bounce evolution. The weakness of their interac- 
tions implies that copious neutrino emission is the 
most efficient means for the hot nascent neutron 
star to cool. Neutrinos diffuse outward through 
the dense core on a time scale of seconds and 
eventually escape with about 99% of the gravi- 
tational energy released during collapse. Some 
of this is redeposited in neutrino captures be- 
hind the stalled shock. But whether or not this 
neutrino heating is a proximate cause of shock 
revival, any definitive account of the explosion 
mechanism must include careful handling of the 
neutrinos, for errors in their overwhelming con- 
tribution to the energy budget reduce confidence 
that an energetically subdominant detail like the 
explosion has been adequately resolved. 

What sort of computation is needed to follow 
the neutrinos' evolution? Deep inside the newly- 
born neutron star, the neutrinos and the fluid 
are tightly coupled (nearly in equilibrium); but 
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as neutrinos are transported from inside the neu- 
tron star, they go from a nearly isotropic dif- 
fusive regime to strongly forward-peaked free- 
streaming. Heating behind the shock occurs pre- 
cisely in this transition region, and modeling this 
process accurately ideally involves tracking both 
the energy and angle dependence of the neutrino 
distribution functions at every point in space. 

2. POST-BOUNCE SIMULATIONS 

A full treatment of this time-dependent six- 
dimensional neutrino radiation hydrodynamics 
problem is a major challenge, and the computa- 
tional resources necessary to solve it are still just 
over the horizon. Nevertheless, much has been 
learned over the years through the simulation of 
model systems of reduced dimensionality. 

Progress on the general trend of increasing to- 
tal dimensionality (space plus momentum space) 
handled by simulations in recent years is sketched 
in Table 1. It should of course be kept in mind 
that this does not represent every insight rele- 
vant to the explosion mechanism or other aspects 
of supernova phenomenology, obtained via sim- 
ulation or otherwise. Nevertheless, the literature 
bears witness that increases in the dimensionality 
of sophisticated simulations are intertwined with 
important advances in the field. (For a previous 
version of this table that includes entries for some 
earlier studies, see Ref. [T].) 

That spherical symmetry would be an incom- 
plete description was suggested both by obser- 
vations of supernova SN1987A and by profiles of 
fluid variables obtained from spherically symmet- 
ric simulations. Working mostly with mock-ups 
of multidimensional instabilities within spherical 
symmetry, one group suggested that a 'doubly- 
diffusive instability' within the neutron star — 
which could boost neutrino luminosities — was 
more important (even crucial) to an explosion 
than convection immediately behind the shock, 
which enhances the efficiency of neutrino heating 
2 3 4 5. Other studies, sacrificing the energy de- 
pendence of neutrino distributions in exchange for 
genuine multidimensional fluid dynamics, found 
the effects of post-shock convection to be more 
important [7J and even robust [6127] . Subsequent 



work could not confirm the presence of doubly- 
diffusive instabilities [9 .37 38 39 , and the opera- 
tion of more conventional convection within the 
neutron star appears to be either suppressed by 
neutrino transport [10] or unimportant to explo- 
sion dynamics [28 39 . And echoing earlier pes- 
simism [5], later work retaining neutrino energy 
dependence also deflated enthusiasm for the ca- 
pacity of post-shock convection to robustly bring 
about neutrino- heated explosions [11128] . 

The nagging qualitative difference between spa- 
tially multidimensional simulations with different 
neutrino transport approximations motivated in- 
terest in the possible importance of more com- 
plete neutrino transport: might the retention of 
both the energy and angle dependence of the 
neutrino distributions improve the chances of ex- 
plosion, as preliminary 'snapshot' studies sug- 
gested 40 22 ? Of necessity, the first such sim- 
ulations were performed in spherical symmetry, 
which nevertheless represented an advance to a 
total dimensionality of 3 (see Table 1). Results 
from three different groups are in accord: Spher- 
ically symmetric models of iron core collapse do 
not explode, even with solid neutrino transport 
[12121123] and general relativity [24126] . Recently, 
however, it has been shown that the more modest 
O/Ne/Mg cores of the lightest stars to undergo 
core collapse (8-10 M & ) may explode in spherical 
symmetry |14|15|16] , 

In an extension to two space dimensions, one of 
these groups deployed their spherically symmet- 
ric energy- and angle-dependent neutrino trans- 
port capability [13] along separate radial rays, 
with partial coupling between rays [29|31j . Ini- 
tial results — from axisymmetric simulations with 
a restricted angular domain — were negative with 
regards to explosions (in spite of the salutary ef- 
fects of convection, and also rotation) [28.31], ap- 
parently supporting the results of Ref. [TT] . An 
explosion was seen in one simulation [29 4 1131] 
in which certain terms in the neutrino trans- 
port equation corresponding to Doppler shifts 
and angular aberration due to fluid motion were 
dropped; this simulation also yielded a neutron 
star mass and nucleosynthetic consequences in 
better agreement with observations than the ex- 
ploding simulations of the 1990s [617] . arguably 
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Table 1 

Selected neutrino radiation hydrodynamics milestones in stellar collapse simulations studying the long- 
term fate of the shock. The 'Yes' entries in the 'Explosion' column are all marked with an asterisk as a 
reminder that questions about the simulations — described in the main text — have prevented a consensus 
about the explosion mechanism. 'Total dimensions' is the average of 'Fluid space dimensions' and l v space 
dimensions,' added to l v momentum space dimensions.' The abbreviation 'N' stands for 'Newtonian,' 
while 'PN' — for 'Post-Newtonian' — stands for some attempt at inclusion of general relativistic effects, and 
'GR' denotes full relativity. A space dimensionality in quotes — like '1.5' — denotes an attempt at modeling 
higher dimensional effects within the context of a lower dimensional simulation. For the fluid, this is a 
mixing-length prescription in the neutron star ('NS') or the heating region ('HR') behind the stalled 
shock. For neutrino transport, it indicates one of two approaches: multidimensional diffusion in regions 
with strong radiation/fluid coupling, matched with a spherically symmetric 'light bulb' approximation in 
weakly coupled regions ('thick/thin'); or the (mostly) independent application of a spherically symmetric 
formalism/algorithm to separate spatial angle bins ('ray-by-ray'). 



Group 


Year 


Explosion 


Total 


Fluid space 


v space 


v momentum 








dimensions 


dimensions 


dimensions 


space dimensions 


Lawrence 


1989 


Yes* 


'2.25' 


'1.5' NS+HR 


1 


1 


Livermore |2|3|4] 








(GR) 




(GR) 


Lawrence 


1992 


Yes* 


2 


2 HR 


2 





Livermore [5] 








(N) 




(N) 


Los Alamos 


1993 


Yes* 


'1.75' 


2 


'1.5' 





EI 








(N) 


thick / thin 


(PN) 


Arizona 


1994 


Yes* 


'1.75' 


2 


'1.5' 





m 








(N) 


ray-by-ray 


(N) 


Florida 


1994 


No 


'2.25' 


'1.5' NS 


1 


1 


Atlantic [El] 








(GR) 




(O(v/o)) 


Oak Ridge 


1996 


No 


'2.5' 


2 


1 


1 


lio|n| 








(N) 




(0(v/c)) 


G arching 


2000 


No, Yes* 


3 


1 


1 


2 


[12113(14115(16] 




(8.8 M ) 




(N) 




(0(v/c)) 


Oak Ridge 


2000 


No 


3 


1 


1 


2 


[17118(19120(21] 








(N) 




(<D(v/c)) 


Arizona 


2002 


No 


3 


1 


1 


2 


22 23j 








(N) 




(0(v/c)) 


Oak Ridge 


2000 


No 


3 


1 


1 


2 


(17118119124125126) 








(GR) 




(GR) 


Los Alamos 


2002 


Yes* 


'2.5' 


3 


'2' 





[6|27| 








(N) 


thick / thin 


(PN) 


Garching 


2002 


No, Yes* 


'3.75' 


2 


'1.5' 


2 (0(v/c), 


13 15 28 29 30131132) 




(11.2 M Q ) 




(PN) 


ray-by-ray 


PN) 


Arizona 


2005 


Yes* 


3 


2 


2 


1 (0(1), 


33 3415] 








(N) 




N) 


Florida 


2006 


Yes* 


'2.75' 


2 


'1.5' 


1 


Atlantic [SB] 








(N) 




(0(v/c)) 
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because of more accurate neutrino transport in 
the case of both observables. The continuing les- 
son is that getting the details of the neutrino 
transport right makes a difference. 

In addition to accurate neutrino transport, low- 
mode (I = 1,2) instabilities that can develop only 
in simulations allowing the full range of polar an- 
gles may ultimately make decisive differences. 

One example is an explosion, obtained with the 
'ray-by-ray' code described above [30115132] . of 
one of the lowest mass stars (11.2 Mq) to have 
an iron core: an explosion was seen when the 
full 180° in polar angle was used, but not when 
the grid covered only 90°. Explosions of 11 Mq 
and 15 Mq progenitors have also been obtained 
by another group on a 180° grid using another 
ray-by-ray code that maintains only the neutrino 
energy dependence (flux- limited diffusion) [36 
These achievements were presaged by an earlier 
study that demonstrated the tendency for con- 
vective cells to merge to the lowest order allowed 
by the spatial domain [42], and especially by the 
more recent discovery of a new standing accre- 
tion shock instability 43J. It is not yet entirely 
clear whether the standing accretion shock in- 
stability simply results from a standing acoustic 
wave in the cavity between the neutron star and 
the shock (44] , or is instead an advective/acoustic 
cycle in the same region [45)46147] ; but in any case 
a suspicion is becoming widely established (e.g. 
48 34 32 36 ) that a shock instability is operative 
in the full complexity of the supernova environ- 
ment that generates asymmetry independently 
of convection [43)49)50] , These shock-instability- 
induced global asymmetries seem to be sufficient 
to account for observed asphericities that have of- 
ten been attributed to rotation and/or magnetic 
fields [5TI35] . 

In another example, this shock instability is ac- 
knowledged but claimed to only set the stage for 
another low-mode instability that is more directly 
related to the explosion, namely, I = 1 (and, later 
and subdominantly, £ = 2) g— mode oscillations 
of the nascent neutron star [34 35 . Excited by 

2 In addition to the shock instability to be mentioned 
shortly, these authors present evidence that the use of a 
nuclear network rather than a representative heavy nu- 
cleus contributes materially to a successful explosion. 



accretion streams funneled onto the neutron star 
by the flow structures resulting from the shock 
instability, the claim is that these oscillations of 
the core radiate acoustic waves that steepen into 
shocks which then deposit sufficient energy at 
larger radii to drive the explosion. This effect 
has so far only been seen by one group, with a 
code whose neutrino transport is two-dimensional 
in space but that retains only the energy depen- 
dence of the neutrino distributions (flux-limited 
diffusion) Jf| It would not have been seen in sim- 
ulations that either exploded by the neutrino- 
driven mechanism before this acoustic mechanism 
had time to develop, or that were stopped 'pre- 
maturely' when it became clear that a neutrino- 
driven explosion was not forthcoming. On the 
other hand, at least one other group !36] has now 
run simulations out to the requisite times with 
a ray-by-ray flux-limited diffusion code, and the 
acoustic mechanism has not been observed. It 
might be claimed that this latter code's use of 
spherical coordinates artificially pins the core to 
the origin and prevents the necessary oscillations 
[35] . However, an order-of-magnitude mismatch 
between the frequencies of the g— mode core oscil- 
lations and the shock instability that would drive 
it is a bit puzzling. Confirmation by other codes 
will be needed to eliminate the possibility of the 
acoustic mechanism being merely an unphysical 
simulation artifact. 

Surely every 'Yes' entry in the explosion col- 
umn of Table 1 has been hailed in its time as 'the 
answer' (at least by some!), and as a community 
we cannot help hoping once again that these re- 
cent developments mark the turning of a corner; 
but important work remains to verify if this is the 
case. At least two groups are pursuing neutrino 
transport in two and three space dimensions that 
retain the full energy and angle dependence of the 
neutrino phase space [52)53)1)54)55] . 

3. SUPERNOVA NEUTRINO SIGNALS 

Whether it be by playing the spoiler, acting as 
agent of explosion, or setting the conditions nec- 
essary to the operation of some other mechanism, 

3 The neutrino transport also does not account for most 
velocity-dependent terms or inelastic neutrino scattering. 
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Figure 1. Crudely estimated neutrino luminosities 
(thick) and characteristic energies (thin). Note the 
change in time and luminosity (but not energy) scales 
between the left and right panels. The left panel is a 
close-up of infall and bounce at t — 0. 

neutrinos have direct or indirect roles in the many 
explosion scenarios discussed over the years — and 
detailed simulation of collapse and the second or 
so after bounce will allow a neutrino signal de- 
tected from a Galactic supernova to serve as an 
unrivaled diagnostic tool. 

But apart from its potential dynamic and diag- 
nostic roles in connection with the explosion, su- 
pernova neutrino emission is of intrinsic interest 
as an energetically dominant feature of stellar col- 
lapse. Figure [TJ obtained from a simple and self- 
contained model [56] . serves as a recognizable if 
very rough caricature of neutrino light curves pro- 
duced by detailed models [31126123157] . Apparent 
are a gradual increase in v e luminosity during in- 
fall, which then drops with neutrino trapping; a 
burst dominated by v e just after core bounce; and 
a cooling phase characterized by emission of v and 
v of all flavors over many seconds. The lower lu- 
minosities and average energies of v e and v e dur- 
ing the cooling phase result from their participa- 
tion in charged-current interactions unavailable 
to the other species at these energies, with v e dif- 
fering from v e because j/ e -capturing neutrons are 
more abundant than i> e -capturing protons. 

At present, detailed simulations focusing on 
the explosion assume massless neutrinos, relying 
on the standard expectation that the effective 
mass from neutrino forward scattering off elec- 



trons suppresses flavor mixing in the high-density 
region where the explosion is launched. Before 
reaching terrestrial detectors, however, neutri- 
nos will pass through flavor-changing MSW res- 
onances at lower densities in the stellar enve- 
lope, and possibly also experience flavor-changing 
effects while traversing the earth. Hence col- 
lapse simulations provide neutrino spectra as in- 
put to studies that take these effects into ac- 
count in determining possible signals in terrestrial 
neutrino detectors. To the extent that generic 
energy- and time-dependent features of the sig- 
nals can be motivated by the large-scale sim- 
ulations, such studies provide a possible basis 
for learning — from neutrino detection of a fu- 
ture Galactic supernova — about remaining un- 
knowns surrounding neutrino flavor mixing. For 
instance, observation of the v e burst could rule 
out sin 2 $i3 > 10~ 3 and a 'normal' mass hierarchy 
[55] , If sin 2 6*i3 > 10~ 3 it may also be possible to 
distinguish between 'normal' and 'inverted' mass 
hierarchies on the basis of whether it is the v e or 
v e signal that is affected by the passage of the su- 
pernova shock wave through the Am 2 3 resonance 
in the supernova envelope [59I6Q] , or by the neu- 
trinos' traversal through the earth [5T]. Passage 
through the earth could also affect both the v e 
and v e channels even if sin 2 #13 < 10 -5 [61j . 

4. NEW EFFECTS AT SMALL Am 2 ? 

Again, the conventional wisdom is that the 
MSW effect suppresses flavor mixing until the 
resonances corresponding to the solar and at- 
mospheric neutrino mass-squared differences are 
reached far out in the supernova envelope; how- 
ever, this standard expectation may in fact 
be wrong. If the off-diagonal contributions 
of neutrino-neutrino forward scattering domi- 
nate, both neutrinos and antineutrinos can mix 
maximally over a significant range of neutrino 
energy — and this could occur deep in the super- 
nova environment, even with the known 'small' 
mass-squared differences [62]. In fact a couple 
of types of collective neutrino mixing behavior 
of potential relevance to supernovae have been 
identified in both analytic and numerical studies 
63 64 65 66J, though their ultimate impact is not 
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yet clear. It may well be that, as a result of non- 
linear effects associated with neutrino-neutrino 
forward scattering, flavor mixing will have to be 
piled on top of the challenge of high dimensional- 
ity that supernova modelers already face in sim- 
ulating neutrino transport. 
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